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The fractal properties of models of randomly placed n-dimensional spheres (n=l,2,3) are studied 
using standard techniques for calculating fractal dimensions in empirical data (the box counting and 
Minkowski-sausage techniques). Using analytical and numerical calculations it is shown that in the 
regime of low volume fraction occupied by the spheres, apparent fractal behavior is observed for 
a range of scales between physically relevant cut-offs. The width of this range, typically spanning 
between one and two orders of magnitude, is in very good agreement with the typical range observed 
in experimental measurements of fractals. The dimensions are not universal and depend on density. 
These observations are applicable to spatial, temporal and spectral random structures. Polydisper- 
sivity in sphere radii and impenetrability of the spheres (resulting in short range correlations) are 
^ also introduced and are found to have little effect on the scaling properties. We thus propose that 

' apparent fractal behavior observed experimentally over a limited range may often have its origin in 

^ , underlying randomness. 



> 

■ I. INTRODUCTION 

(N 
^ . 

, In recent years the study of fractal structures has been an active field of research both theoretically and experimen- 
tally In theory, a variety of algorithms and dynamical models which produce fractal sets have been introduced. 

Typically, in these models one can define an asymptotic limit in which the set exhibits fractal behavior on an arbi- 
trarily broad range of length scales. One can then approach this limit by a process of gradual refinements of the set, 
which may involve either an increase in system size or decrease in the minimal object size. In the case of empirical 
fractals observed experimentally the situation is different. For these fractals the range over which they obey a scaling 
law is restricted by inherent upper and lower cutoffs. In most experimental situations this range may be quite small, 
namely not more than one or two orders of magnitude [|lO|. Nevertheless, even in these cases the fractal analysis 
(-H ■ condenses data into useful relations between different quantities and often provides useful insig ht d l]. 
Q I Motivated by the yet largely inexplicable abundance of experimentally observed fractals, we consider in this paper 
O ' the apparent fractal properties of systems which are governed by uniformly random distributions. The choice of 
^ random systems is justified by the abundance of randomness in nature, and by the fact that uniform randomness is a 
convenient limit, on top of which correlations can be introduced as perturbations. Although a purely random system 
cannot be fully scale invariant, it may, as we show below, display apparent fractality over a limited range. The width 
and the cut-off values of this range seems to be in good agreement with the typical width and typical cut-offs observed 
in experimental measurements (between one and two decades) , unlike the case with models which are inherently scale 
free. 

To illustrate these ideas we consider a model in which n-dimensional (n = 1,2,3) spheres of diameter d are randomly 
distributed in an n-dimensional space in the regime of low volume fraction occupied by the spheres. In three dimensions 
(3D) our definitions coincide with ordinary spheres, while for n = 2 (2D) we consider disks, and for n = 1 (ID) rods 
of length d. In the basic model the positions of the spheres are uncorrelated and they are thus allowed to overlap. 
We then extend the model to the case where there is a distribution of sphere radii and examine the effect of this 
distribution on the fractal properties. We also examine a version of the model in which spheres are not allowed to 
overlap, thus introducing short range correlations between sphere positions. This class of models may approximately 
describe spatial distribution of objects such as craters on the moon, droplets in a cloud and adsorbates on a substrate. 
In particular, the one dimensional model may describe the level distribution in energy spectra of quantum systems 
and the zeroset of random temporal signals. Therefore, such models may be at the root of empirical observations of 
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fractals in experiments and data analyses, dealing with processes governed mainly by randomness |I0 . As fractality 
is usually revealed by applying various resolution analyses, we also address the question of whether our observations 
are method-dependent. 

Two of the most commonly employed resolution analysis methods are the box-counting (BC) and Minkowski sausage 
(MS) techniques [|l],fll . The fractal properties of these models are studied here both analytically and numerically, 
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within the BC and MS frameworks. The analytical solution is exact except for the case of impenetrable spheres in 
n > 1 dimensions where some approximations were needed. 

In the box-counting (BC) procedure one covers the embedding space by non-overlapping boxes of linear size r, 
and then counts the number of boxes N(r), which have a non-empty intersection with the (fractal) set. A fractal 
dimension Dbci is declared to prevail at a certain range, if a relation of the type 

iVsc(r) -r-^«^ (1) 

holds, or equivalently, if the slope of the log- log plot 

Dbc = -slope {logr, \og[NBc{r)]} (2) 

is found to be constant over that range. In the MS case one draws an rt-dimensional (n=l,2,3) Minkowski-sphere 
(M-sphere) of radius R around each point in the set under consideration, and calculates how the volume V of the 
union of all spheres changes under a variation of R. The set is considered fractal with a dimension Dms i over a range 
of scales, if 

y(i?) - (3) 

or equivalently, if 

Dms = -slope {logi?, log[y(i?)/i?"]} (4) 

is constant within this range. The BC and MS methods are known to be identical from the mathematical point of 
view , where the limits r, R ^ can be taken. Their equivalence from the physical point of view, under the 
constraints of finite cut-offs, is not obvious. While it is clear from Eqs.(|^)-(|^) that both methods employ resolution 
analysis, in which the number of occupied "resolution units" (boxes or spheres) has to be determined as a function 
of the resolution magnitude (box-length or M-sphere-radius), there are slight differences between the two methods, 
which are due to both geometrical differences and the presence of cut-offs. 

Using a resolution analysis one should be aware of the existence of physical cut-offs. This is especially important 
in the structures considered below which are not fractal in the rigorous mathematical sense. We will show that the 
log-log plots of the functions N{r) and V{R) display linear behavior between these physical cut-offs. The slope in 
this range can be interpreted as an FD. The existence of this scaling behavior between physically relevant cutoffs is 
a central motivation for the study of the random models presented below. 

The paper is organized as follows: in Section || we consider the basic model of penetrable spheres where all the 
spheres are of equal size, and calculate the BC function. From the analysis of this function we obtain the nontrivial 
linear range and extract the apparent FD. The generality of this model is then demonstrated using an information 



theory argument. In Section III we generalize the analysis to the case where the spheres are not equally sized but 



exhibit various size distributions and examine the effect on the apparent fractal behavior. Certain correlations are 



then introduced in Section IV where the case of impenetrable spheres in considered and the conclusions are given in 
Section 0. 



II. THE MODEL OF RANDOMLY DISTRIBUTED SPHERES 



A. The ID Model 

In this model M rods of length d <C 1 are randomly placed on the unit interval (Fig.Q) such that the 
positions of the rod centers arc chosen from a uniform random distribution. The rods are mutually penetrable, 
namely overlaps are allowed. The BC function N{r), will now be derived. This function gives the number of boxes, 
for given box size r, which have a non-empty intersection with the set. For a large enough number of rods (of the 
order of 100 in the present case), the deviations from the expectation value {N{r)) are negligible and edge effects are 
unimportant. Let p denote the probability that a box of size r intersects a rod of length d. Then for a total of 
boxes, 

{N{r)) = P. (5) 
r 

Following Refs. [ p"3[p^ ] define the probability qi that after random placement of the first rod, a given box remains 
unintersected. Neglecting edge effects, this implies that the center of the rod must be at least a distance of d/2 away 
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from either edge of the box (Fig.|l|(c)). A total length of r + 2{d/2) is, therefore, unavailable for placement without 
intersection. For a uniformly random distribution it follows that: 



qi = l-{r + d). (6) 

The next rods are placed independently, which means that after placement of M rods, the probability that the box 
is still unintersected is 

<Z-(l-(r + d))*^. (7) 
Finally, the probability of at least one intersection is p = 1 — (7, and therefore: 



1 



r 



(iV(r)) = 5 l-{l-{r + d)f . (8) 



This is the one dimensional BC function for randomly adsorbed, mutually penetrable rods | [lO| . Numerical simulations 
of the model along with the theoretical prediction of Eq. (||), are shown in Fig.^. The excellent agreement is evident. 
The parameters d and M are independent, and are limited only by the restriction r/ < 1, where rj = Md is the coverage 
of the line. We find that for large M (here M > 100) changing d and Af while keeping 77 fixed, merely translates 
rigidly the BC function in the log-log plane. 

In FigJI, the BC function is shown in a range between cut-offs, together with a linear regression. This is the range 
which is commonly measured in experimental work. The special significance of these cut-offs is discussed below. 
Notice the nearly linear appearance, extending over close to 1.5 decades. 



B. Extension to 2D and 3D 



In the two dimensional model one places disks of diameter d on the unit square, so that for a given box of 
area to remain unintersected, no disk center may fall within the area shown in Fig.|l](d). Thus, a total area 
of +4:{j'K {d/2)'^) + 4:{rd/2) is excluded for placement of the first disk center. Therefore in this case the probability 
of the box to remain empty is: 



qi = l- l^r 

The next disks are placed independently, leading to: 



2 ' -Trd^ + 2rd 
4 



mr)) = \ 



1 - ( 1 - ( + 2rd+ ^Trd^ 



A/ 



(9) 



(10) 



Similarly, for 3D one considers independent placement of spheres, and the excluded volume for sphere centers is that 
formed by convolution of a cube of side r and a sphere of diameter d, which is + 6(r^ d/2) + 12(r jT:{d/2)'^) + 



[| |7r((i/2)^). One obtains then 



91 = 1 



4 6 



and 



{N{r)) = ^ 



1- ( 1- ( r3 + 3r2d+— rd2 + ^d3 
4 6 



M 



(11) 



(12) 
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C. Cut-offs, Range of Linearity, and the Fractal Dimension 

1. Cut-offs 

We will now study the BC function N{r) and examine the possibility of fractal-like behavior. For simplicity we will 
first concentrate on the ID case, but the conclusions apply to 2D and 3D as well. 

For mathematical fractals displaying full scale-invariance, the log-log plot of N{r) vs. r can form a straight line 
with a fractal slope over an unlimited range of scales in the asymptotic limit. For the set we consider, this is clearly 
not the case as in FigJ^: there are "knees" beyond which the slope approaches 1. This is due to the existence of lower 
and upper cut-offs, tq and ri respectively. These cut-offs correspond to relevant physical limits of observation. Here, 
the smallest feature is of size d, and so the finest resolution is of that size. No additional information is obtained by 
reducing r below vq = d , where D approaches the trivial limit of 1 as r —s- 0. To see why this is so, suppose, for 
convenience, that the center of each rod is located at a point connecting two adjacent boxes: then halving the box 
doubles the number of intersected boxes if r < d, with the result that D must approach the limit of 1 as r ^ 0. When 
the centers of the rods are located arbitrarily, the lower cut-off will not be sharply located at d. Nevertheless, 

ro^d (13) 

is a good estimate for it. As for the upper cut-off ri, it is the size beyond which practically all boxes intersect at least 
one rod, where again D 1. This happens when the boxes are larger than the average gap between rods: 

ri^l/M-d, (14) 

which is therefore an approximate upper cut-off. We thus have the approximate range ro < r < ri, where the 
measurement is properly tuned to measure inherent scaling behavior, if it exists. For scaling behavior to be observed, 
there must be a minimal range of apparent linearity (Eq.(^)). The size of such a range and the extent of linearity 
displayed by the BC function in it are considered below. 

2. Range of Linearity 

The standard experimental procedure is to apply a linear regression analysis on the log-log presentation of the 
scaling range. The linear regression line is constructed to go through the inflexion point (r^, N{ri)) of the log-log plot 
of Eq.(||). The dependence of the range of linearity A on the coeflicient of determination, TZ^, (measuring the quality 
of the linear description) is then explored. 

The range of linearity is approximately given by A = 2(log(ri) — log(ro)). Applying a linear regression analysis on 
the log-log plot of the BC function, we evaluated the slopes and actual ranges of linearity under these constraints, 
with different values of TZ'^ imposed. Typical results are shown in Fig.^. For instance, about two decades of linear 
behavior can be observed for a required value of TZ^ of below 0.97. Examples of experimentally observed fractal 
objects exhibiting several orders of magnitude of linearity are rare; the vast majority of reported fractality spans 1-2 
orders of magnitude ||^-|^, [lo| , |l5| . It is important to emphasize that, in fact, we are mimicing in our calculations the 
common practice of searching and reporting FD's in empirical data. 

The convention for the smallest meaningful scaling range is one decade ||l6|| . Given this, a simple argument yields 
the maximum allowed coverage: By using the estimates above for the cut-offs, ro = d and ri — I /M — d = d(l — 77) /rj, 
one observes that a range A — log(ri) — log(ro) = log[(l — r])/r]] of at least one decade, requires that: 77 < 1/11. It 
follows that the apparent fractality in our model is restricted to rj < 0.1. The range between the cut-offs grows as the 
coverage is decreased, as observed in Fig.^. 

In addition to the width of the range between the cutoffs, the quality of the linear fit within this range, measured 
by TZ^ should also be considered. One can limit the range of linearity by imposing a lower bound on TZ^: obviously, 
the range decreases as TZ'^ increases (Fig.^). Also note from this Figure, that for a given range, the quality of linearity 
grows as the coverage is increased, or as the slope of the BC function between the cut-offs (i.e., the FD) approaches 
1 (FigH). This happens because of the smooth merging with the slope beyond the cut-offs, which is trivially 1. We 
thus conclude that the two cut-offs limit the width of the linear range for high coverage while the TZ^ criterion limits 
it for low coverage. As a result, the range of scales in which we observe apparent fractality is typically between one 
and two orders of magnitude. 
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3. Fractal Dimension 



The apparent FD shown in Fig.||, for an imposed values of TZ^ — 0.995, rises monotonically from as more rods 
are added. This is an important aspect of the model: it does not predict a universal (specific) FD, but the whole 
allowable range of FD values. The regression results are further compared in Fig|| to an analytical equation, obtained 
by calculating the logarithmic derivative of N(r) at the estimated middle point Ve = ^/rori: 



D 



(1) _ 

BC — 



d{\og[N{r)]) 
d(log[l/r]) 



1 - V?Ki - v)- 



M 



M-1 



1-1- 



M 



M ■ 



(15) 



We use an estimate for the middle point {r^) rather than the exact result (r^), since cannot be given analytically. 
The almost symmetric S-shape of the log-log plot in the scaling region assures that is a good estimate for r^. As 
seen in FigH, the FD predicted by Eq.([l5|) is an accurate lower bound to the regression result. 

By using (1 + x/N)^ — > as iV ^ oo, Eq.(15) may be simplified in the "thermodynamic limit", while keeping 77 
finite. One then obtains: 



D 



(1) 

BC 



exp 



v/^(l - '7) 



?7+ V'7(l-??) 



(16) 



Notably, the expression for -D^^ depends on 77 alone. In the limit of small 77, one can further simplify Eq.(16) and 
obtain 



D 



(1) 

BC 



V 



1/2 



77 < 1. 



(17) 



Expressions for -D^^ and -D^^ may be derived from Eqs. 
discussion of their cut-offs and range of linearity are deferred to Sec.|IIA2|, where a more general model is treated. 



0^,(fr3) for the 2D and 3D cases. These, as well as a 



D. Absence of Correlations 



Since fractal objects typically exhibit some correlations, one might wonder whether the finding of an apparent FD 
in our model, is also due to some hidden correlations within a certain window of resolution. For instance, it might 
seem plausible that the finite extent of the rods introduces a correlation, for if a point on the line belongs to a certain 
rod, then a point at a distance x < c? is likely to belong to the same rod. However, correlations at a scale x < d below 
the lower cut-off are, rightly so, typically not measured. Furthermore, as we show next, the correlation exponent for 
X > d vanishes. This will prove that the emergence of an FD in the penetrable rods case is not due to the presence 
of correlations, but entirely due to the scaling displayed by the underlying uniform distribution. 

A correlation exponent, is associated with a given set if a correlation function c{x) has the following power-law 
form over a sufficiently large range: 

c{x) - x"". (18) 

An exact expression for this function may be derived for penetrable rods. Let 0{xq) be the local density at the point 
xq on the hne. That is, 9{xo) = 1 if xq belongs to a rod, and zero otherwise. The correlation function is defined as: 

c{x)^{9{xo)e{xo+x)) (19) 

where (• • •) denotes either an ensemble average (at fixed xq), or an average over xq. Assuming ergodicity, we choose 
to calculate the latter. Only pairs of points (xq, xq -I- x) such that both 9{xo) and 9{xo + x) are 1 contribute to the 
average. We thus require the simultaneous probability: 

P(x, xo) = Pr mxo +x)^l)n (0(xo) = 1)] . (20) 

The correlation function is then: 
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c{x) — / dxo 6{xo) 9{xq + x) P{x,xo), (21) 
Jo 

where one can set 0{xo) = 6{xo + a;) = 1. Neglecting edge effects, P{x,xo) does not depend on xq, since we are 
considering a statistically homogeneous medium (or stationary process). In what follows, therefore, Xo is considered 
to be any convenient reference point, though none of the results depends on its location. With this general choice of 
xq, one obtains from Eq.(^l]): 

c{x) = P^„ (x) = Pr [{e{xo +x) = l)n {e{xo) = 1)] . (22) 

To evaluate Pxg (x) , consider a modification of the argument that led to Eq. (Q) : pick two "test-points" at random on 
a line of length L and denote their respective positions xq and xq + x. Now consider randomly placing rods of length 
d on the line. The events which are complementary to both points {xq,xq + x) being occupied, are (1) that at least 
one is unoccupied (with probability Qi), and (2) that both are unoccupied ((52(2^)). Clearly Qi already accounts for 
Q2{x), so that: 

P,„{x) = l-[2Qi-Q2{x)]. (23) 

The evaluation of Qi is a repetition of the argument leading to Eq.(^, with a vanishing box size (r = 0). Thus: 

Qi - (1 - d)"*' ^ e-". (24) 

when M — > 00. Q2ix) requires that no rod center is placed within a distance smaller than d/2 to either test-point 
(see Fig.|l|(e,f)). When x > d, this clearly excludes a total length of 2d from possible placement of rod centers. When 
X < d, a, length of only x + d is excluded, due to overlap of the two inner segments. In all other respects the argument 
leading to Eq.(|^) (with r = 0) is again repeated, so that: 

Q2ix) = (1 - (d + six)) f ^ e-''(''+^(^», (25) 
where = 1/M is the average distance between rod centers, and 

<-^={t\:i'd. (26) 

Combining Eqs.(p^p6|), one finally obtains for the correlation function: 

, , r 1 - 2e-" + ■, ^ > d 

'^^^^ ~ \ l- 2e-'^ + e-'^e-P'' : x < d . ^ 

Thus for X > d, c{x) is constant, i.e., the correlation exponent vanishes (a result obtained numerically by Meakin and 
Deutch [p^). However, the non-trivial range of the similarity dimension starts at d, and hence cannot possibly be due 
to correlations. These exist, not surprisingly, for x < d, but even there they decay exponentially (with a characteristic 
length of 1/M), and not as a power-law, as required for fractality measured by the scaling of correlation functions. 

This completes the demonstration that the non-trivial result we obtained for the similarity dimension of penetrable 
rods is not due to correlations. The elements of this demonstration are applicable also to the more general model 
presented in the next section and are not repeated there: the absence of correlations holds for all penetrable cases 
treated in this study. 



E. Generality of the Model 



In this section arguments are presented showing that the model of randomly placed, mutually penetrable spheres 
introduced in Sec.^ is very general in the sense that it is a prototype of a much larger family random processes. 
Being minimalistic in assumptions, the model is in fact a generic one for random processes. An information-theoretical 
approach is now employed to show this, using the fact that the only assumption entering the model is the knowledge 
of a mean quantity. An important feature of the information-theoretical approach is that it extends the model from 
spatial to temporal random sequences, and even to energy-level distributions JTsf . Thus, in an adsorption process 
the relevant mean quantity is the average density of adsorbates and one seeks the distribution of nearest-neighbor 
distances. An equivalent situation is that in which one knows the average period of a time signal and is interested 
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in the distribution of intervals between successive crossings of the time axis (the "zeroset"). A finite width d is then 
the signal width at the crossings. Yet another case is that of spectral-level distribution (where d is the uncertainty 
band width) , for which it is well known that spacings of quantum energy levels in classically non-chaotic systems are 
exponentially distributed |l9[| . Indeed, it will be shown that an exponential distribution of the intervals between the 
centers of the spheres is an inherent characteristic of our model. It shall then become clear that spatial, temporal and 
spectral random processes fall into the same class. The spatial and temporal cases differ only in that a time-process 
must be ordered on the time axis, while the positions of "events" on a spatial axis are not necessarily ordered. There 
are, however, iV! possible arrangements in the spatial case, obtained by interchanging the labels of events on the 
spatial axis. The time-process may be considered as the one ordered set out of all these permutations. Relabeling is 
all that is required to map a spatial process to one in time, and clearly this does not affect the statistics of positions 
or intervals. Furthermore, it is well-known |^ that the FD of the actual time signal can be deduced from that of 
its zeroset. Relying on the equivalence just discussed, we choose in what follows for concreteness to work on the 
spatial process of placement, but one should keep in mind that the discussion applies just as well to time-processes 
and energy- level spacings. 



1. Information Theory Argument 

Consider now the derivation of the distribution satisfying the assumption of knowledge of the mean. In the present 
case this is the known average placement density 

M 

{x) = - (28) 

(or equivalently, this may represent the average period of a time signal.) The arguments presented here for ID are again 
easily generalized to 2D and 3D. Following a standard information theory argument pO| ], the resulting distribution 
P{x) of distances x between neighboring adsorbate centers is obtained by maximizing the missing information, S, 



5 = - / P{x) \nP{x)dx, (29) 
Jo 

where the constraint of knowing (x) can be written as: 

/"OO 

{x) = / xP{x)dx. (30) 
To this one must add the normalization constraint 

/•oo 

/ dx P{x) = 1. (31) 
Jo 

Using Lagrange multipliers A and the maximization of S can then be written in terms of a functional F as: 

/■oo 

F[P{x)]^- dx{P{x) In P{x) + XP{x)+fixP{x)}, (32) 
Jo 

whose variational derivative is: 

/>oo 

6F^- dx{\nP+l + X + fix}SP. (33) 
Jo 

The arbitrariness of SP then requires the vanishing of the term in curly brackets, or: 

P{x) = e-(i+^)e^^^. (34) 
Inserting this into the constraint equations, Eqs.(p0[),(|3l|) yields the distribution of nearest-neighbor distances 



P(x) = ^e^Tiy. (35) 

The appearance of an exponential distribution is not surprising: it is the Maxwell distribution when the average 
energy is given as the constraint in the canonical ensemble; or it shows up as the distribution of time-intervals 
between successive radioactive decays, where the mean life-time acts as the known constraint. It remains to be shown, 
returning to the adsorption-language, that given this exponential distribution of intervals between the centers of the 
rods, the adsorbate positions are uniformly distributed. 
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2. The Position Distribution 



It is a standard exercise in probability theory to show that a uniform distribution of positions leads to an exponential 
distribution of intervals (see e.g. Ref. or recall the argument for the time of flight of a particle undergoing random 
collisions with a given mean free path). It is now shown that the opposite holds as well, namely that the exponential 
distribution of intervals, derived above from an information-theory argument, leads to a uniform distribution of 



positions. The argument that led to the derivation of the exponential distribution of intervals, Eq.(35), assumed that 
only the mean distance between points is known, and that this is the only parameter of relevance. Therefore, it was in 
fact implicitly assumed that successive placements are independent (for otherwise additional constraints should have 
appeared, reflecting the dependence of the distribution of intervals on the number of previous placements). Let P{x) 
denote the probability density of finding a point between x and x + dx after a single placement. Given that there is a 
point at xq, consider the conditional probability density g{xo + a;|a;o) of finding the nearest neighboring point between 
xq + x,xq + x -\- dx. This can be expressed as: 

g{xo + x\xo) ^ -j^e^i^ P{xo + x)dx, (36) 

where eyip{—x/{x))/{x} is the probability density of finding a gap of length x. But since this exponential probability 
depends only on the (non-negative) distance between neighboring points, it is clear that nothing prevents interchanging 
the roles of xq and xo + 2:, i.e., it must hold that: 

g{xo + x\xo) = g{xo\xo + x), (37) 



or, explicitly: 



1 1 X 

— e P{xq + x)dx = — e (-) P{xo)dxo. (38) 

(x) (x) 



One is at liberty to choose dx = dx^, so that P{xq) = P{xq A- x), which holds for every x. Therefore P{x) is constant^ 
i.e., the positions are uniformly distributed. This result followed from the exponential distribution which was derived 
under the minimalistic assumption of knowledge of the mean of a relevant property. Drawing on the generality of 
this derivation, we conclude that a uniform random distribution of adsorbate centers is a generic model of random 
processes in space and time, where one only assumes knowledge of the mean. A uniform distribution of adsorbate 
centers is, however, exactly what was assumed in the adsorption model in Sec.||, for which apparent fractality was 
detected. Fractality may therefore be expected for any other random system which can be characterized by its mean. 

One may further employ the above information-theory formalism in order to derive the distributions appropriate 
to knowledge of higher moments, if correlations are present in the system. 



III. POLYDISPERSED MUTUALLY PENETRABLE SPHERES 



A. The Model 



The basic model introduced in Sec^is now generalized, by allowing polydispersivity in radii. That is, we consider a 
model of randomly placed n-dimensional spheres with a distribution of radii P{a). The radii are assumed to be chosen 
independently from P{a). In ID, this may, e.g., represent a random spectrum with levels exhibiting a distribution 
of lifetimes. In 2D, such systems may approximately describe, for example, the formation of metal clusters on metal 
surfaces |^-^. In 3D one might consider the distribution of atmospheric or intergalactic dust aggregates. Thus 
the polydispersed case represents a very wide class of systems, whose possible scaling properties and apparent FD 
are quantities of interest. Both the BC and MS functions and FDs will be calculated, first generally, and then for a 
number of specific but broadly used radii distributions. 

1. Minkowski Analysis 

The Minkowski function and dimension corresponding to the model of polydispersed spheres introduced above, is 
calculated next. 
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a. The Minkowski Function 

Let all lengths be normalized to the total linear extent of the surface L. For a given probability distribution of 
radii P{a), consider the determination of the MS function for M randomly placed spheres with a specific realization 

of radii out of P{a,), a = (ai,a2, ...jOm)- The MS function for this realization, V:^^\r), is the volume of the imion 
of spheres of radius R centered at all points in the set under consideration. This amounts to increasing the intrinsic 
radius of a sphere Oi to Oj + R, and then calculating the volume of the union. Let the total volume of the embedding 
space be Vt- It is convenient to work with the normalized volume = Va"'\R)/Vi. The calculation of a^\R) 

is very similar to the calculation of the BC function. A point is randomly chosen in the embedding space, and one 
calculates the probability, qm, that after placement of M spheres with radii a, the chosen point is not included in 
the volume of any of the spheres. The probability qi for this to happen after random placement of the first sphere is 
proportional to the volume remaining after subtracting the volume of this sphere. 



gi = l- 



7„(ai+it!)" 
Vt 



where 



2 : n = 1 

7„ = .( TT : n = 2 



^ : n = 3 



(39) 



(40) 



is a geometrical factor associated with the volume of an n-dimensional sphere. The next placements are independent, 
so that 



M 



= n 1 - 



i=l 



Vt 



(41) 



Now, I — qm is the probability of finding the chosen point in the set, which on the other hand is equal to the normalized 
volume of the set. Thus: 



M . 

ai")(i?) =i-ni 

i=l ^ 



M 



(42) 



where p = M/Vt is the number density. 

Next, average over all possible realizations of radii. Let 



(/(a))^ jp{a)f{a) 



da 



(43) 



denote the average of any function of the radius. Then the expectation value of the normalized volume, when each 
realization of radii a is weighted by its probability Hi^i P{'^i) dcii, is given by: 

M 



aW(i?) = (ai")(i?)) = 1 - / []da, P(a,) 1 - 



Plnjaj + Py 
M 



y..P(., (i _ ^2^)] " = , _ (i _ ,2.ii^y- , 



In the limit of large M one finally obtains: 

Q(")(i?) = l-e-''^''<(''+'^)">, 

and, explicitly for each dimension: 

a(^) {R) = l- exp [-2p{R + (a))] 

a^'^^R) = 1 - exp [-np{R^ + 2R{a) + (a^))] 

a^^\R) = 1 - exp -^piR^ + 3i?^(a) + 3R{a'^) + (a^)) 



(44) 



(45) 



(46) 
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b. The Fractal Dimension 

To find the FD according to Eq.(^, we evaluate the logarithmic derivative of Q;*-"-'(i?), which is found to be: 



dlog 

d\og{l/R) 



l-exp(-p7„((a + i?)«)) 



(47) 



In analogy with the discussion in Sec. [I C , one may expect the cut-offs to be found approximately at 

1 



compared to ro = d and ri = p 



Ro^{a) and i?i = _ (48) 
— d for BC. The difference from BC by a factor of 2 is due to the appearance 



of radii as opposed to box- lengths (which are equivalent to diameters). As in Sec. [I C 2, we define the FD to be the 
slope at the estimated middle-point _Re = VRqRi of the scaling range: 



D 



(") _ 

MS ~ 



dlog {a^''\R)lR'') 
d\og{l/R) 



(49) 



which yields: 



D 



(1) _ 

MS ~ 



1/2 ,1/2 



exp {rii + yjriiil - r/i)^ 



/ , 1/2 jl/2\ 

exp \rii + u^' (j)^ J 



D 



(2) _ 
MS ~ 



3/4 ,1/2 
^2 <P2 ■ 



O 1/2 , 



, 3/4 ,1/2 , 1/2 , 
exp 772+t^2 <f'2 +^2 (P2 



1 



^(3) 
^MS 



= 3 



(AV3 ^ _ (2^)2)1/% 2i.3^/^ 03 + 3z.3^/' 0^ 



exp + (A2/3 - (2^^)2)^/2 + 03 + ^.1/2 03/2) 

where for convenience a number of parameters are defined as follows: 

Vn= Pin (a") 

(a)" 



(50) 



P 
A 



?7n 
1 

2 
3 

2^73 



l/n _ l/„ 
In '^n 

{a') (a) 



V3- 



3 

73 (a3) 



2\^3/2 



(51) 



The parameter 77,1 is simply the coverage p^ ] ; the other parameters measure various moments of the radius distribution. 

The expressions for D^y^ can be simplified somewhat for the case of a constant radius, whence (a") — a", and 
consequently 7y„ = v^- Note that the FD contains useful information about the first n moments of the radius- 
distribution. In the next section the BC function of the same model of polydispersed radii is solved for, in order to 
compare the MS and BC methods. 



2. Box-counting Analysis 

a. Tfie BC Function 

Once again, take all lengths to be normalized to the total linear extent of the surface L. Repeating the averaging 
arguments used to arrive at Eq. (^^ , combined with the derivations of Eqs. (p|) , ([To|) , (p^ , leads to the following result 
for the BC (density-)function in the polydispersed case: 
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M 



fe=0 



Here r is the normalized box-length, p ~ M/Vt is the number-density, and 

Pkn — 



(52) 



In 


: k = 





1 


: fc = 


n 


2n 


: fc = 


n — 




: fc = 


1, n 



(53) 



is introduced for convenience. Taking the limit of large M, one finds: 



7V("Vr) = — 



1 — exp 



XI f^kn r'' {a 



n—k\ 



(54) 



Although at first sight this r esult m ay seem to differ substantially from that for the MS function a'"-'(i?) (Eq.^5|), the 
two will be compared in Sec. Ill B 2 and shown to be quite similar. 

b. The Fractal Dimension 

Proceeding with the calculation of the FD, evaluation of the logarithmic derivative yields: 



dlog(iV(")(r)) p(ELifc/?fcn^'(«""^)) 
dlog(l/r) exp (p J2k=o f^kn r'=(a"-'=)) - 1 ' 



(55) 



As usual, the FD is defined as the slope at the estimated middle-point re = y^roTi between the lower and upper 
cut-offs, To — 2(a) and ri = — 2(a): 



D 



(„) _ dlog(iV(")(r)) 



BC 



dlog(l/r) 



(56) 



with the result: 



^(1) 

^BC 



= 1 



\/f?i(l-»7i) 



^1/2 ,1/2 



exp {rji + V??i(l - - 1 exp (rji + ■0^^^ 



(2) 
BC 



^BC 



,1/2 



= 3 - 



((A')^/3;x-(3M)^)^/%CrV'3 + |C3^/^Vr 



exp (,3 + ((A')^/3 M - (3m)^) + ^3 + K3'' ^3'') 
where we defined for convenience: 

6n = 



Cn 

A' 



- 1 



(57) 




2/ 2XN3/2 



((37r)>2» 
73(0^ 



(58) 



Here (5„ is a geometrical factor associated with the volume of an n-dimensional cube (or box). Note the similarity of 
Cnjtpn, A' to Vn,4'n,^ respectively of the previous section. 

The BC result for the one-dimensional case is identical to the MS result. This is due to the fact that in ID both a 
box and a spher e reduc e to a line segment. For n = 2, 3 the geometrical factors are different (7„ for the MS and (5„ 
for BC). In Sec. [II A3 these issues will be considered in detail. 
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3. The Distribution Functions 



In this section the effect of polydispersivity is examined exphcitly, by assuming various functional forms for the 
radii distributions. 

To assess the influence of polydispersivity we considered four types of common continuous distributions of radii, 
a well as a simple discrete bi-modal distribution. These are compared for reference with the case of monodispersed 
radii treated for BC in Sec.|^. The distributions considered are: 

1 -<°-<°>)^ 

Pn(o,) = /— , e 2p Normal 
Pu{a) = + Uniform 

PE{a) = (iye""/^''^ Exponential [26] (59) 

Ps{a) = ^^b-i-^/'a^e-/" Schulz || 

PB{a) = ( l'^;" = {"j 7 Bi - modal 

^ ' \ Pr[a — [a) + p b\ = 1 — p 

where b is in all cases a measure of the width of the distribution. A comparison of the distributions for two combinations 
of (a) and b values is shown in FigJq, where they can be seen to differ significantly. In order to understand the difference 

(n) (n) 

between these distributions in the present context, note that -D^(A and Dj^/g are determined by the first n moments 
of the distributions. Therefore it is useful to summarize the differences as follows: let y„ — (a")/(a)" and z = b/ (a). 
Then: 

y2 = 1 + z'^ 2/3 = 1 + Normal 
^2 = 1 + 12 y3 = l + x Uniform 

2/2 = 2 2/3 = 6 Exponential (60) 

2/2 = 1 + 2 2/3 = 1 + 3z + 2z2 Schulz 

y„ = p(l - (1 -p)z)" + (1 -p)(l +pz)" Bi-modal 

However, the effect of assuming these different radii-distributions on the MS and BC functions, and on the respective 
FDs is marginal in spite of the differences among them, as shown next. 



B. Results 



1. Range of Linearity and the Effect of Dimensionality 

In order to meaningfully compare the results in different dimensions, it is most convenient to fix the average distance 
(x) = p~^/^ between sphere centers. Note that this implies different coverages, since from Eg. (pT|) : ?]„ = 7„((a) / (a;))". 
In particular, since typically (a) ^ (a;), 7/3 will be much smaller than r/i for the same average distance. With this choice, 
a general estimate for the range of linearity, independent of the dimension, may be found. Suppose (x) = 10*^ (a); 
using the values for the cut-offs, i?o = (a) and Ri = ip^^/" — (a), one has: 



A„ < logi?i - logi?o = log {\{x)/{a) - 1) = log (ilO'^ - l) - A: - 0.3 for fc > 1. As discussed in Sec. |ll C % the range 
of linearity is limited both by the distance between cut-offs, which tends to increase when the coverage decreases and 
the quality of the linear regression, measured by 7?.^ which tend to improve as the coverage increases. These trends 
are independent of the embedding space dimension and therefore in all dimensions we observe apparent fractality 
within a range of 1-2 orders of magnitude. 



2. Comparison of Box- Counting and Minkowski Sausage Results 

It is comforting to find that, by and large, the BC and MS methods of resolution analysis yield very similar results. 
This is shown in Fig.^ where the two methods are compared for the monodispersed case (quantitatively similar 
differences result for poly-dispersivity). When properly normalized (see caption), the methods are identical in ID, 
and differ slightly in 2D and 3D. We shall see below that the FD values do not differ by more than 0.05 for 2D and 0.1 
for 3D either, with MS giving the consistently lower values (see Figs. ^ ^ below). We attribute the small variance 
between the two methods to the finite size of the basic building blocks, and to the differences in the geometrical 
factors 7ra and (5„, which determine the details of lowering the resolution of observation. As the difference is so small, 
in the next section the entire discussion is held in terms of BC. 
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3. Effect of Polydispersivity on the BC Function 



The BC function is displayed as a function of yardstick size in Figs.^,^for the distributions considered above, for a 
representative density of spheres, corresponding to an average intersphere distance (x) of 20 sphere radii (see captions 
for details). The average radius (a) and width 6 were set equal for all distributions. The value of 6 = 0.9(a) was 
chosen so as to reflect a very broad distribution. Values of b larger than (a) are not permissible since this would 
lead to negative radii. An observation of significant importance is that the qualitative (and to a large extent also the 
quantitative) nature of the results is unaffected by changing h. The BC function for the ID case depends only the 
average radius, as is clear from Eqs.(^),(^, and hence no difference is observed between distributions with the same 
(a) in the ID case. For 2D {a?) enters, which depends, in turn, on the specific distribution (Eq^p9|)). For 3D also {a?) 
enters, so that a stronger dependence on the distributions results. However, as seen in Figs. the differences in BC 
function for different distributions set in only at small r values, close to the lower cut-off, and are always bounded. 

The two most important questions in the context of apparent fractality relate to the slopes and the ranges of 
linearity. The former is dealt with in detail in the next subsection. As for the range of linearity of the scaling region, 
it appears that this is slightly decreased when polydispersivity is compared to monodispersivity. 

4- Effect of Polydispersivity on Fractal Dimension 

The FD is displayed as a function of coverage in Figs. |l0| , |ri| for the distributions of Eq.(^9|). The average radius (a) 
and width b were set equal for all distributions, with b = 0.9(a) chosen again so as to reflect the unfavorable case of 
a very broad distribution (see captions for details). Changing among distributions is seen to have only a minor effect 
on the FD, for the entire coverage range. A somewhat stronger effect is observed for 3D than for 2D. For ID, as 
), the FD depends only on the average sphere radius (i.e., rod length), which was taken identical 
for all distributions, so that they necessarily all produce the same FD. However, the second and third moments (a^), 
(a'^), which determine Z?^^) and D'^'^^ through Eqs.( p0|),(|57| ), do depend on the particular distribution, and must affect, 
therefore, the FD. Comparing the moments of Eq.(|60|) reveals why the FD is so robust with respect to change of 
radius-distribution or its parameters: recalling that z < 1 in order to assure positivity of the radii, observe that ?/„ are 
typically close to 1 (the monodispersed case) in all cases, except for the exponential distribution. (For the Bi-modal 
distribution it is straightforward to show that < ?/„ < 1 (recall that < p < 1) , a nd is close to zero only for very 
small p and z.) However, even for the exponential distribution, it is seen in Figs that although the corresponding 
FD is indeed somewhat displaced, it is still very close to that of the other distributions. Thus, it is the combined effect 
of the relatively small dependence of the moments on the underlying radius-distribution in the parameter range of 
interest, together with the even further suppressed sensitivity of the FD to these dependences, which is responsible for 
the robustness of the FD. The insensitivity of the model to polydispersivity is in support of our proposition that the 
random adsorption model is generic: its features are virtually unaffected by (strong) perturbations in this commonly 
encountered way. 

IV. SOLUTION OF THE IMPENETRABLE SPHERES CASE 

A. The Model 

Consider now a different model, which adds correlations on top of the model of equi-sized, mutually penetrable 
spheres by imposing impenetrability on a system of n-dimensional spheres at equilibrium . Impenetrability creates 
a negative correlation in sphere positions. This model is fully solvable for n = 1, and approximately solvable with high 
accuracy for n = 2,3. It represents an important class of processes with correlations, such as models of hard-sphere 
liquids, energy-level repulsion in quantum systems which are classically chaotic pj| , Langmuir-type adsorption, etc. 
As demonstrated below, the correlation due to impenetrability merely modifies the apparent fractal character already 
induced by the random nature of the problem. We limit ourselves to a BC analysis of this problem. 

B. Derivation for ID Case 

Consider first the ID case of impenetrable rods: BC function is calculated by employing a result derived using 
thermodynamic arguments by Helfand et al. |^^, and using a statistical argument by Torquato et al. They 



seen in Eqs.(|q),(57 
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show that the probabihty of finding a cavity of length I containing no rod centers, in a system of a large number of 
impenetrable rods of length d each, is: 



^(1) 



(1 - ?7i)exp 



1 



(61) 



for / > d. In performing BC, for a box of length r to remain empty, no rod-center may fall within d/2 from either 
side of the box, so that I ~ r + d. g*^^^ may be rewritten, in the limit of a large number of rods {M ^ oo at fixed 
coverage), as: 



1-^1 



M 



Having p ~ 1 — q and using Eq.(||) the expected number of intersected boxes is: 



{N{r)) 



l-(l-77i) 1- 



(62) 



(63) 



1. Fractal Dimension 

As in the penetrable rods case, one can now use Eq.(||) (with the slope calculated at r = r^) to calculate a lower 
bound for the FD. We obtain: 



,1 'sJt^m 



This can again be simplified for large M: 



D = l- . (65) 



C. Derivation for 2D and 3D Cases 



So far all the results were based on exact calculations. We now consider an approximate solution for 2D and 3D 
impenetrable spheres. Full analytical solutions are at present impossible: exact results for the probability of finding 
a cavity containing no disk or ball centers after their placement, as employed in Eq. (|6l|) a bove, are not available 
because the n-particle probability densities are not exactly known. Nevertheless, Refs. [30 3^ provide some accurate 
approximations for the probability of finding a two- or three-dimensional spherical cavity of radius / in an equilibrium 
system of hard spheres of radius a. 



1. BC Function 



In the context of BC, one in fact requires the probability of finding a cavity with the shape of the convolution of 
a box and sphere (see Sec.|l| and Figj^(d)). This will be undertaken in a future study; at present we will settle for 
an approximation of the cavity by spheres. The important quantity to conserve in this approximation, is the cavity 
volume, since this is what actually enters the probabilistic argument at the root of the calculation of the BC function. 
If one simply takes the spherical cavity radius as the geometric mean of the sphere radius plus half box length, and 
sphere radius plus half box diagonal. 
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the real cavity volume (see Eqs.(p|jTl|)) is overestimated by no more than 10% in 2D, and (except for a sharp maximum 
of 80% for 0.2 < r/a < 5) by no more than 20% in 3D. This will suffice for the present purpose of an approximate 
treatment of the 2D and 3D impenetrable spheres problem. Now, the results of Ref. for the empty-cavity 
probability in n dimensions can be expressed conveniently as follows: 

g(")(r) = (1 - r;„) exp (^^^5^/„(r)) , (67) 

where the subscript n on the coverage 77 serves to remind that the same coverage defines different combinations of 
sphere sizes and densities for different dimensions (see Eq.(pTl) in the monodispersed case). The functions fn depend 
on the particular approximation used, but have the general form 

/«« = E«?'(^»)<; x„ = ^. (68) 

The a^"'' are, for n = 2 (impenetrable disks), using the scaled-particle theory of Reiss, Frisch and Lebowitz ||3^ : 

=277-1 ; ^ = -Ar] ; = 4, (69) 
whereas for n = 3, using the Carnahan-Starling result [Q, Ref. ||3^ finds: 

af^ = -^(977^ - 77; -f 2) ; ^ = 127?^ ; a^^^ = -677(8 + t]) ; = 8(1 + 77). (70) 
Next, the BC functions are given, as usual, by: 

(N^"Hr)) = ^{l~q^''^), (71) 
The apparent FDs predicted by these results are discussed next. 



2. Fractal Dimension 



For the 2D and 3D cases, the FD is obtained by evaluating the logarithmic derivative at ri^^ = [2a - 2a)] 

Using Eq.(pl[), the estimated middle-point of the scaling range can be rewritten as: 



1/2 



= 2ft.„(77) ; h„{r]) 

a 



1 / \ 1/*^ 

1 / 7ri 



1/2 



2 \rin 

The hard-sphere radius a cancels out and the FD for tt, = 2, 3 is found to be given by: 

(„) d(log[jV(")(r)]) , 1 Vn ^< rtoh ( ^^n^iW 

^ = d(log[l/r]) 1^=- ^ " - ^JT^^^hMiV^i^hM + 1) + 1] X 



(72) 



exp 



2-^'4"^('7«)[(l + hM){^ + V^KiVnW^ - (1 - Vn) 



(73) 



It should be noted that Z?*^"^ are functions of alone, which indeed, in contrast to the penetrable spheres case, is 
exactly the volume fraction of space occupied by the impenetrable spheres. 

D. Results 

1. Effect of Impenetrability on the BC Function and on the Range of Linearity 

Fig.|l^ shows the plot of Eq.(||) for the penetrable rods case , together with the expression for {N{r)) in the ID 
impenetrable rods case, Eq.p3. The behavior is qualitatively similar in both cases, and virtually indistinguishable for 
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low coverages. Fig.|l2| thus demonstrates the primary role of pure randomness in the appearance of fractality, even in 
the presence of correlations, at least at coverages below 10%. 

The impenetrability case results for 2D and 3D are shown in Figs.|l^ and ^ respectively, along with those for 
penetrable spheres. As in the ID case, the effect of impenetrability is virtually unnoticcablc in terms of the FD. Once 
more, it appears that the scaling content of this highly non-trivial system is already contained to a very large extent 
in its penetrable counterpart. 

2. Effect of Impenetrability on Fractal Dimension 

The result for ID is shown in Fig.||. As expected, the FD rises to 1 faster than in the penetrable rods case: when no 
overlap is allowed, the line is filled up at lower coverages. The important observation is the large range of coverages 
for which the penetrable and impenetrable rods FDs nearly overlap. For all practical purposes, therefore, the FD is 
very close in these two cases, which differ significantly in the extent of correlations present in the respective systems: 
the fractal content of the impenetrable rods geometry is already contained to a large extent in the associated random 
penetrable rods case. 

V. CONCLUSIONS 

We have shown in this paper that randomness in its most elementary forms, generates apparent fractal structures 
over 1-2 decades, between a lower cut-off - the elementary building block - and an upper cut-off which is approx- 
imately the average distance between building blocks. We adopted an empirical-like approach to the calculation of 
fractal dimension: we deliberately considered sets bound within finite cut-offs, which display scaling behavior over 
a physically relevant range. We believe this approach to be both useful and necessary, if a direct contact between 
theory and experiment is to be achieved. This led us to consider several simple, but widely applicable models of 
random phenomena, with and without correlations, obtaining as one of the main results, an analytical solution for 
the apparent fractal dimension of models of randomness. The models studied are convenient starting points for other, 
more elaborate ones of random phenomena. It is argued that the models we have studied are, in fact prototypical of a 
large class of spatial, temporal and even spectral random phenomena. The methods introduced here should be useful 
in the study of further cases of random phenomena, with other types of disorder and correlations. To conclude, we 
would like to point out that apparent fractal behavior is expected to occur for a more general class of distributions. 
Generally, the logiV(r) vs. log(l/r) plot includes lines of slope D = n (n = 1, 2, 3) beyond the upper and lower 
cut-offs which are connected by an interval of slope D < n, which depends on the specific distribution. 

Acknowledgements 

Wc would like to thank P. Bak, R.B. Gerber, C. Jayaprakash, D. Mukamel, J. Sethna and G. Shinar for very helpful 
discussions. D.H. and D.A. are members of the Fritz Haber Research Center for Molecular Dynamics. D.A. is also a 
member of the Farkas Center for Light Energy Conversion. 



16 



URL: 
URL: 



http : // www.fiz .huj i . ac . il/ st aff _ Acc / faculty /biham 



*** E-mail: davidOgranite.fh.huji. ac.il 
[1] B. B. Mandelbrot, The Fractal Geometry of Nature (Freeman, San Francisco, 1982). 

[2] Fractals in Physics, Essays in Honour of B.B. Mandelbrot, edited by J. Feder and A. Aharony (North Holland, Amsterdam, 
1990). 

[3] Fractals in Science, edited by A. Bunde, S. Havlin (Springer, Berlin, 1994). 

[4] On Growth and Form, No. 100 in NATO ASI Ser. E, edited by H. E. Stanley, N. Ostrowsky (Martinus Nijhoff, Dordrecht, 
1986). 

[5] The Fractal Approach to Heterogeneous Chemistry: Surfaces, Colloids, Polymers, edited by D. Avnir (John Wiley & Sons 

Ltd., Chichester, 1992). 
[6] H. Takayasu, Fractals in the Physical Sciences (J. Wiley & Sons, Chichester, 1990). 

[7] Scaling Phenomena in Disordered Systems, No. 133 in NATO ASI Ser. B, edited by R. Pynn, A. Shjeltrop (Plenum Press, 
N.Y., 1985). 

[8] P.-G. de Gennes, Scaling Concepts in Polymer Phystcs (Cornell University Press, Ithaca, 1979). 

[9] A.-L. Barabasi and H.E. Stanley, Fractal Concepts in Surface Growth (Cambridge University Press, Cambridge, 1995). 
[10] D.A. Hamburger, O. Biham, D. Avnir, preprint, 1995. 

[11] K. Falconer, Fractal Geometry: Mathematical Foundations and Applications (Wiley, Chichester, 1990). 

[12] For a line of length L, normalize all lengths to L. 

[13] H.L. Weissberg, J. Appl. Phys. 34, 2636 (1963). 

[14] S. Torquato, G. Stell, J. Chem. Phys. 79, 1505 (1983). 

[15] In Ref. [0, we include a histogram, based on all the experimental FD measurements presented in Ref. [^, and showing 
the number of plots as a function of the number of decades of the linear range. The conclusion is that by and large, most 
experimental measurements of fractal dimensions are based on data that extends between one and two decades. 

[16] P. Pfeifer, M. Obert, in The Fractal Approach to Heterogeneous Chemistry: Surfaces, Colloids, Polymers, edited by D. 
Avnir (John Wiley & Sons Ltd., Chichester, 1992), pp. 38-39. 

[17] P. Meakin, J.M. Deutch, J. Chem. Phys. 83, 4086 (1985). 

[18] L.S. Cederbaum, E. Haller, P. Pfeifer, Phys. Rev. A 31, 1869 (1985). 

[19] M.V. Berry and M. Tabor, Proc. R. Soc. A 356, 375 (1977). 

[20] F.M. Reza, An Introduction to Information Theory (McGraw-Hill, New York, 1961). 

[21] K.A. Brownlee, Statistical Theory and Methodology (John Wiley & Son, New York, 1965), p. 169. 

[22] H. Roder, E. Hahn, H. Brune, J.-P. Bucher, and K. Kern, Nature 366, 141 (1993). 

[23] H. Brune, C. Romainczyk, H. Roder, and K. Kern, Nature 369, 469 (1994). 

[24] G. Rosenfeld, A.F. Becker, B. Poelsema, L.K. Verheij, and G. Comsa, Phys. Rev. Lett. 69, 917 (1992). 

[25] ri„ is an appropriate variable in the context of sequential addition of n- dimensional spheres, such as in an adsorption 

experiment. When one is a priori faced with the final result of overlapping spheres, such as in random spectra or the 

distribution of craters on an asteroid (M. J. S. Belton et ah, Science, 265, 1543 (1994)), it is more convenient to work with 

the net coverage, q'"'(0) = 1 — exp (?;„), of Eq.(^). 
[26] A nearly exponential distribution of relatively compact (i.e., circular) clusters was found in embedded-atom simulations of 

Ag on Pt(lll) (P. Blandin, C. Massobrio, and P. Ballone, Phys. Rev. Lett., 72, 3072 (1994)). 
[27] The Schulz distribution is a continuous version of the discrete Gamma distribution. Its standard formula is -Ps(a) = 

/ \ m + l 

T(^Ji+i) (l^J a™e-(™+i)"/<''> (G. V. Schulz, Z. Phys. Chem. B 43, 25 (1939)), where m > -1 determines the 

peakedness of the distribution. A general expression for the nth moment is: (a") — b" YYi-i (^T^ + • parameter b we 
introduced is a measure of the width of the distribution: b = (a) — a* , where a* = (a) ^^-^ is the position of the maximum 
of Ps (a). 

[28] Equilibrium and non-equilibrium distributions such as random sequential adsorption (RSA) are known to be different at 

the same coverage; see B. Widom, J. Chem. Phys., 44, 3888 (1966). 
[29] E. Helfand, H.L. Frisch, J.L. Lebowitz, J. Chem. Phys. 34, 1037 (1961). 
[30] S. Torquato, B. Lu, J. Rubinstein, Phys. Rev. A 41, 2059 (1990). 
[31] S. Torquato, S.B. Lee, Physica A 167, 361 (1990). 

[32] H. Reiss, H.L. Frisch, and J.L. Lebowitz, J. Chem. Phys. 31, 369 (1959). 
[33] N.F. Carnahan and K.E. Starling, J. Chem. Phys. 51, 635 (1969). 



17 



FIG. 1. Typical realizations of random placements of rods of size d (gray with center dot) in the penetrable (a) and impene- 
trable (b) cases. Also shown is a division into "boxes" (vertical lines). The figure also represents the zeroset of a temporal series 
of signals of width d, or a spectrum of randomly positioned energy levels with uncertainty-width d. (c) Illustration of counting 
procedure used in our probabilistic arguments; For a box of length r to remain unintersect ed, n o rod center may approach its 
ends closer than a distance of d/2. Hence a total length of r -|- d must remain free. For Sec. [IB: (d) Excluded area in the case 
of placement of disks in 2D. The excluded area consists of the box of side length r, four quarter circles of diame ter d , and four 
rectangles r ■ d. This is just the convolution of the box and a circle of diameter d. (e) Counting procedure in Sec, II D: Distance 
between the "test points" x is larger than d: for both points xo and xo + x to be empty, a total length of 2d must be excluded, 
(f) If a:: < d, a rod falling in between the two points may overlap both, so a length of only x + d is excluded. 



FIG. 2. Comparison of simulation results (circles) to the theoretical prediction of Eq.(|8|) (solid line) for the number of 
intersected boxes as a function of their size in the ID penetrable rods case. The coverage is rj — 0.1 and the rod length is 
d/L = 10"^ The cut-offs are manifested as the two knees in the graph. The lower bound ro is seen to be indeed located at 
r — d. The upper bound ri is at r = d/77 — d, also conforming with the prediction in the text. The agreement between theory 
and simulations is excellent over the entire range. Inset: Same with 77 = 0.01 and d/L = lO"'^. Note the increase in the range 
of linearity. 



FIG. 3. The simulation results (circles) for the number of intersected boxes N{r) vs. r in the experimentally relevant range 
are shown along with a linear regression fit for coverage t] = 0.1 (obtained for d = 10~^, Nd ~ 10'*). This is the experimentally 
relevant range which is typically used to obtain the FD. 



FIG. 4. The range of linearity. A, of Eq.(g), in a linear regression analysis. The range of linearity decreases as higher quality 
regression is required (see text). The results presented are valid in all dimensions, but it should be remembered that the same 
coverage corresponds to different inter-particle distances in different dimensions. 



FIG. 5. Apparent fractality (FD) as computed by linear regression with a relatively high (0.995) coefficient of determination, 
in the case of rods. The prediction of the analytical equations, Eqs.(^) and ( |65| ) serve as accurate lower bounds. The 
prediction for the FD of the penetrable and impenetrable rods cases is seen to differ only marginally, indicating that the 
dominant contribution to the FD comes from the penetrable (i.e., totally random) rods case. The lowest coverage shown 
corresponds to the lowest molecular densities observed in nature: lOA/cm"^ in intergalactic space. 



FIG. 6. Probability distributions P{a) used for polydispersivity in radii of penetrable spheres model: Normal, uniform, 
exponential, Schulz, bi-modal. All distributions have the same average (a) and width b, defined in the text. Shown here is the 
wide case, (a) = 1, b — 0.9. 



FIG. 7. Comparison of box-counting (log[A'^(r)] vs log(r)) and Minkowski (log[a(_R)] vs log(_R)) functions for monodispersed, 
mutually penetrable spheres in n = 1,2,3 dimensions. The linear density is fixed at pi — 0.05/a, where a = 1 is the 
sphere-radius. A meaningful comparison is achieved by plotting the Minkowski function normalized to the sphere volume 
{'jnR"), with R — r/2. In ID BC and MS are identical. Differences do develop, albeit for ruler sizes (r, R) beyond the cut-offs, 
for n = 2 and n = 3, which depend on (a^) and (a^), (a^) respectively, and have different geometrical factors due to the use of 
boxes and M-spheres. Note that N{r) and a{R) are normalized to the total number of boxes and the total volume respectively. 



FIG. 8. BC functions for poly-dispersed, mutually penetrable spheres in n = 2 dimensions. As in Fig.^, the linear density is 
pi = 0.05/(a) (note that p has dimensions of inverse area). The distribution parameters are (a) = 1, b = 0.9. As seen clearly in 
the inset, the exponential distribution lies highest, followed by Schulz, normal, bi-modal, uniform, and finally the monodispersed 
case (no distribution). Since the lower cut-off is expected at log 1 — 0, it appears that the effect of polydispersivity is to 
decrease the range of linearity somewhat. 



FIG. 9. Same as Fig.g, but in n = 3 dimensions (p has dimensions of inverse volume). The asymmetry between the lower 
and upper cut-offs is apparent. 
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FIG. 10. Fractal dimension, obtained as slope of the BC and Minkowski functions at the estimated middle points r^, Re, for 
polydispersed, mutually penetrable spheres in 2D, as a function of the coverage, logj^Q ri2. The minimal coverage is as in Fig.^. 
Displayed here are the broad-distribution results: (a) = 1, b — 0.9(a). Small differences are observed among the distributions. 
The order is opposite to that in Fig.^, with the Minkowski dimension consistently somewhat smaller. Also shown is the FD 
for impenetrable disks, which as expected, is slightly larger. The lowest coverage shown corresponds to the lowest molecular 
densities observed in nature: lOA/cm^ in intergalactic space. 



FIG. 11. Fractal dimension, obtained as slope of the BC and Minkowski function in 3D. The details are the same as in Fig.[lO| 
and the results are qualitatively the same as in 2D. 

FIG. 12. Comparison of box-counting predictions in penetrable and impenetrable rods cases. The results for penetrable 
(Eq.(^)) and impenetrable rods (Eq.(|63|)) virtually coincide for rj < 10~^. For tj = 0.1 a barely noticeable difference develops. 
d = 10~® in both cases. 
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